Aging and Crossovers in Phase-Separating Fluid Mixtures 
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We use state-of-the-art molecular dynamics simulations to study hydrodynamic effects on aging 
during kinetics of phase separation in a fluid mixture. The domain growth law shows a crossover 
from a diffusive regime to a viscous hydrodynamic regime. There is a corresponding crossover in 
the autocorrelation function from a power-law behavior to an exponential decay. While the former 
is consistent with theories for diffusive domain growth, the latter results as a consequence of faster 
advective transport in fluids for which an analytical justification has been provided. 
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I. INTRODUCTION 



There has been intense recent interest in understand- 
ing the coarsening dynamics of phase-separating mix- 
tures [H4T3j. For fluid systems, attention has focused 
on single-time quantities, e.g., the correlation function 
and structure fac tor, or the growth law of the domain 
scale £{t) @, [13, H EEj]- An equally important and 
deeply related aspect of coarsening systems is the na- 
ture of the aging properties [l6| which are encoded in 
two-time quantities. This issue has been intensively 
studied in the context of disordered systems and struc- 
tural glasses, in phase-separating systems without hy- 
drodynamic modes, such as solid mixtures, but not in 
segregating fluids. In part, this is due to the scarcity 
of reliable numerical results for these systems. In this 
letter, we present the first molecular dynamics (MD) 
results for aging in coarsening binary fluid. Before dis- 
cussing our model and results, it is useful to give a brief 
overview of the main concepts. 

Following a quench inside the miscibility gap, a ho- 
mogeneous binary mixture (A+B) separates into A-rich 
and B-rich phases. The evolution of the system from the 
randomly-mixed phase to the segregated state is a com- 
plex process. Domains rich in A- and B-particles form 
and grow nonlinear ly with time [H, Q3, Gil, GJ] • This 
coarsening is a self-similar process which is reflected in 
the scaling behavior of various physical quantities, as 
the two-point equal-time correlation function G(r, t) — 
(0(fi,i)0(f2,t)) - {(f>(ri,t)){(f>(r 2 ,t)) (where the order 
parameter </> is the difference between the A and B con- 
centrations and r = \fi — |)- This quantity can be 
expressed in the scaling form G(r, t) — G{r/£), G{z) be- 
ing a time-independent master function. Typically, £(t) 
grows in a power-law manner as £(t) ~ t a , where the 



growth exponent a depends upon the transport mecha- 
nism. For diffusive domain growth |19j a = 1/3, which 
is referred to as the Lifshitz-Slyozov law. This behav- 
ior characterizes segregating solid mixtures. However 
in fluids, hydrodynamic transport mechanisms, impor- 
tant at large length scales, lead to much faster growth. 
In spatial dimension d — 3, for a percolating domain 
morphol ogy , advective transport results in three growth 
regimes [2C 




e(t) ^e dv = (D V ) 1/2 , 

Idv <£l(t) <&£ in = if /(fry), 
£(t) » t in . 



(1) 



In Eq. ([T]), D is the diffusion constant, -q the viscosity, p 
is the density, 7 the surface tension, idv is the crossover 
length from diffusive to viscous regime, and £j n th e one 
from viscous to inertial hydrodynamic regime [l5| ■ 

For diffusive phase-separation, also the order- 
parameter auto- correlation function C{t,t w ) = 
{4>(f,t)4>(f,t w )) — ((f)(r,t)){(f)(r,t w )) (where t and the 
waiting time t w < t are two generic instants) takes 
a scaling form C(t,t w ) = C(£/£ w ) where hereafter 
£ = £(t) and £ w = £(t w ). The scaling function behaves 
as C(x) = x~ x for large x = £/£ w , A being the 
Fisher- Huse (FH) exponent [2l| whose value in d = 3, 
according to the results in Ref. [22j], is in the range 
[1.5,1.6]. Despite this good understanding of the 
aging properties in diffusive systems, the behavior 
in segregating fluids remains largely unknown. In 
this letter we start filling the gap by studying the 
properties of the autocorrelation function in a range 
of times encompassing the first two regimes of Eq. 
P in a model of binary fluid. In doing that we 
uncover an interesting structure with a crossover which 
corresponds to the switch from the early diffusive to 
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the intermediate hydrodynamic regime in Eq. ([T]) for 
the growth law. Specifically, we find the two parameter 
scaling form 

C(t, 

iw j ^dv 

) = C(x,y), (2) 

where x — £/£ w , y — idv/iw, and C(x,y) obeys the 
limiting behaviors: 

r( r „\ - J X ~ X (with A ~ 1.5), j/>l; 
°^'W-\ e -(*-i)/r 2/ «i. W 

Thus, in the diffusive regime [l w <C ^du or y 1), 
we obtain a power-law decay of the scaling function as 
for spin systems. To the best of our knowledge, this 
has not been observed in previous studies of fluid phase 
separation. Furthermore, and more importantly, in the 
hydrodynamic regime (£ w ^> £a v or y -C 1), we find 
a crossover to a novel exponential. The latter finding 
is the central result of this letter. Here y ^> 1 and 
y <C 1 should be read as nearly pure diffusive and pure 
hydrodynamic regimes, respectively. 

The rest of the paper is organized as follows. We 
describe the model and methodology in Section II. Re- 
sults are presented in Section III. Finally, Section IV 
concludes the paper with a brief summary. 

II. MODEL AND METHOD 

We consider a periodic box of linear dimension L con- 
taining A and B particles of equal mass m. Particles i 
and j, at distance r apart, interact via the potential 

f u(r) = U(r) - U(r e ) - (r - r c ) (f ^ ; r < r c , 
\ u(r) = 0; r > r c . 

(4) 

Here, the Lennard- Jones pair potential U(r) has the 
form U(r) — Ae a p[(a/r) 12 — (cr/r) 6 ], where a is the in- 
teraction diameter, and e a p is the interaction strength 
between particles of species a and f3 [a,/3 G (A,B)]. 
The potential in Eq. (HI) is cut-off at r — r c , to ensure 
faster computation (23l [2^|. For convenience, we set 
r c = 2.5a. The third term on the right-hand-side of 
Eq. (|4]), after a shifting of the potential by its value at 
r = r c (see the second term), was introduced to make 
both the potential and force continuous at r — r c . The 
overall density p is set to unity so that the fluid is in- 
compressible. For the choice eaa = £bb — 2eab = e, 
we have an Ising-like fully symmetric model. The phase 
behavior and equilibrium properties of this system, that 
exhibits a liquid-liquid phase transition at a critical 
temperature ksTc ~ 1.423e, are well studied _25j. For 
convenience, we set e, ks, m and a to unity below. 



Our choice of a high density ensures that the gas-liquid 
transition is well separated from the liquid-liquid one. 

We employ this model to study the kinetics of phase 
separation via MD simulations [23|,[24[ by quenching ho- 
mogeneous critical (50% A and 50% B particles) config- 
urations, prepared at a very high temperature (T = 10), 
to temperatures below T c . We have used the Verlet ve- 
locity algorithm with integration time step At — O.OIt, 
where the time unit to = (ma 2 /e) 1//2 = 1. The temper- 
ature was controlled via application of a Nose-Hoover 
thermostat [24l |, that is well-known to preserve hydro- 
dynamics. 

While unambiguous confirmation of expected hydro- 
dynamic effects, via MD simulation, in binary fluid 
phase separation was done only recently, by using this 
model 0, Hfij], signature of such fast hydrodynamic 
mechanism was observed in a number of earlier MD 
simulations [27l - l29j . 

For analysis of our results, we mapped the contin- 
uum fluid configurations onto a simple cubic lattice of 
size L 3 0]. There each site i is occupied by a particle 
and we assign a spin value Si = +1 if it is an A-particle 
and Si = — 1 elsewhere. This mapping ensures that the 
pattern dynamics can be studied in a manner analo- 
gous to the spin-1/2 Ising model. The two-point equal- 
time correlation function was calculated as G(r, t) = 
(Si{t)Sj{t)) - (SiityiSjit)), with r = \i - j\ and the 
angular brackets denote an averaging over independent 
runs. For a conserved order-parameter system, G(r,t) 
exhibits damped oscillations around zero. The average 
domain size £(t) is obtained from the first zero crossing 
of G(r, t). The auto-correlation function was calculated 
asC(t,t w ) = {Si(t)Si(t w ))-{Si(t)}(Si(t w )}. All results 
presented subsequently for T = 1.1 = 0.77 T c corre- 
spond to L = 64, whereas those for T = 1.25 = 0.88 T c 
were obtained for L — 48. For statistical quantities, we 
average over 10 independent runs. 



III. RESULTS 

Figure Q] shows snapshots from the evolution of the 
binary fluid. The snapshot at t = corresponds to 
the homogeneously mixed phase immediately after the 
quench to the temperature T = 1.1. The other snap- 
shots show the growth of bicontinuous A-rich and B- 
rich domains - this continues until the system reaches 
equilibrium. 

To obtain a quantitative understanding of the evolu- 
tion, in Fig.[5]we plot t(t) for two different temperatures 
[2| . The inset shows the collapse of data for C(r, t) from 
different times upon scaling the abscissa by £(t). The 
behavior of £(t) is consistent with the expected diffusive 
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FIG. 1: Evolution snapshots obtained after a quench to a 
temperature T = 1.1. The A-particles are marked black, 
B-particles are gray. 
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FIG. 2: Growth £(t) as a function of time for temperatures 
T = 1.1 and T — 1.25, on a log-log scale. The solid and 
dashed lines correspond to diffusive (a = 1/3) and viscous 
hydrodynamic (a — 1) growth laws, respectively. Inset: 
Scaling plot of G(r,t) from three different times for T = 1.1. 



behavior (with a = 1/3) at early times, with a crossover 
to a faster growth later on. As discussed in [26|, the 
second regime is consistent with linear hydrodynamic 
growth at later times. Indeed the slightly off-parallel 
nature (on a log-log plot) of the late-time data from 
the linear behavior is due to non-zero off-sets at the 
crossovers that can be conveniently subtracted off to 
obtain a genuine linear behavior [2g]. We observe that 



the crossover between the two growth regimes is seen 
to occur earlier for lower temperatures, since diffusion 
is enhanced by thermal fluctuations. 
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FIG. 3: Plot of C{t,t w ) vs. x = £(t)/£(t w ) for T = 1.1 
and different values of t w (Main: t w ~ 10, 500, Inset: t w = 
10,20,30,40,50), on a log-log scale. The solid lines denote 
the power-law behavior x~ 3 ^ 2 . 

As stated earlier, our primary objective is to study 
aging phenomena, essentially the behavior of C(t,t w ). 
In Fig.[3j we plot C(t, t w ) as a function of x = £(t)/£(t w ) 
for different t w . It is observed that, for the smaller val- 
ues of t w , presented in the inset, there is a reasonable 
collapse of data, not perfect though, and the master- 
curve obeys the form C(x) ~ x~ x , with A ~ 1.5. This 
is the expected behavior for phase-separation in the dif- 
fusive stage. We must add here the reason for the im- 
perfect collapse of data which will also explain the logic 
behind the choice of a narrow range of t w to demon- 
strate this power-law scaling behavior. The height at 
the beginning of the power-law regime changes with t w , 
the increase of which takes the domain order-parameter 
closer to the equilibrium value. Further, with the in- 
crease of t w , faster decay due to hydrodynamics (see 
below) shows up earlier in x. Due to this latter fact the 
data sets as a whole give the impression of an exponent 
steeper than the expected value. However, from the re- 
sults for t w — 10 in the main frame, consistencey with 
the FH exponent should be appreciated. 

On the other hand, for larger values of t w , e.g. 
t w = 500 in the main frame, a crossover from the FH be- 
havior is observed, where, in the post crossover regime 
the decay of C(t,t w ) is much faster than a power law. 
In order to check for a possible exponential decay at 
large t w , a semi-log plot is presented in Fig. 0|a) . In 
this figure, the values of t w were judicially chosen such 
that for time beyond this, diffusive mechanism is prac- 
tically negligible compared to the influence of hydrody- 
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FIG. 4: (a) Same as Fig[3] but on a semi-log plot and for 
t w = 1000, 1500, 3000. (b) As in (a) but at T = 1.25. 



namics. The smallest value of t w in this figure, viz., 
t w = 1000, was chosen being motivated by a finite-size 
scaling study [26[ that showed a pure viscous hydrody- 
namic behavior for time beyond this. Interestingly, in 
this post-crossover regime a scaling form is recovered, 
as proven by the excellent data collapse. Moreover, the 
behavior of the scaling function C (x) is very nicely con- 
sistent with an exponential decay. A similar behavior, 
with scaling and exponential decay of C(x) is also ob- 
served for T = 1.25 (Fig.HJb)), showing that the whole 
pattern is rather general. Here we note that the rate of 
decay, defined by r, in this hydrodynamic regime should 
depend upon the interfacial tension and the advection 
field (influenced by transport properties like shear vis- 
cosity). This will be justified towards the end of this 
paper. As is clear from parts (a) and (b) of FigJH this 
decay rate is a temperature dependent quantity. 

The deviation from a linear look, observed for T = 1.1 
at very large values of x, occur in a range of time where 
finite-size effects have been observed [26|. Indeed, as it 
can be recognized in Fig. [TJ for such times the size of 
domains is comparable to the system size, the liquid is 
close to equilibrium and thus the configurations change 



very slowly. For this reason, we are unable to access 
the inertial hydrodynamic regime [a — 2/3 in Eq. ([I])] 
which would require significantly larger system sizes and 
is beyond our present computational resources. 

We attribute this exponential decay of the autocor- 
relation function to a fast advective field in the hydro- 
dynamic regime which is explained below. The contin- 
uum dynamical equation [Cahn-Hilliard] , used to study 
kinetics of diffusive phase separation in binary alloys, 
can be modified to [lal 



DV 2 n, 



(5) 



to investigate hydrodynamic effects in fluid phase sepa- 
ration. Here v is the time and space dependent velocity 
field that also has dependence upon shear viscosity as 
well as other transport properties and /i is the chemical 
potential. In the regime where hydrodynamics is the 
dominant mechanism, the right hand side of Eq. ([5} 
(that takes care of diffusion) can be neglected to write 



dt 



+ v- V0 = O. 



Assuming v to be a constant and representing 
fc-space, one obtains the solution 



(6) 
in the 



4>{k,t) = <f)(k,t w ) 



a iv-k{t — t w ) 



Hence for the fc-space autocorrelation C(k,t) 
{(f>(k,t)4>(—k,t w ))], one has 

C(/c, i, t w ) = G(k, t w ^j e* ^ w , 



(7) 



(8) 



where Q{k 1 t w ) — C(k, t w , t w ) is the equal-time structure 
factor for time t,„ . The autocorrelation function then is 



C(t,t w ) = / dkC(k,t,t w ) 



= G(f=v(t-t w ),t w ),(9) 

where G(r,t w ) is the (real space) correlation function 
at time t w . Recalling the scaling property G(r,i) = 

G (j^y) ■ one arrives at 

'v(t - t w y 



C(t : t w ) — G 



£(t w ) 



(10) 



For small values of x, Porod's law [15j states that 
G(x) ~ 1 — ax, where a is a constant. This is consistent 
with an exponential decay for small x. The question 
then is what happens for large values of x. Here note 
that an analytical form of G(x), for conserved order- 
parameter dynamics, still remains a challenging task. 
However, our analysis of the numerical data [see the 
inset of FigJ2] is suggestive of an exponential decay of 
the two-point equal time correlation function for large 
enough values of x. 



5 



IV. CONCLUSION 

We have undertaken the first study of aging dynam- 
ics during phase separation in a binary Lennard- Jones 
fluid. The average domain size grows in a power-law 
fashion, with an exponent a = 1/3 (diffusive) at early 
times and a = 1 (viscous hydrodynamic) at later times. 
This crossover has a remarkable consequence for the 
scaling form of the autocorrelation function C{t, t w ). In 
the diffusive regime, C(t,t w ) shows a power-law decay 
with the Fisher-Huse exponent. However, in the hy- 
drodynamic regime, we observe an exponential decay. 
This is due to the fact that advective hydrodynamic 
flows wash out the memory very rapidly by displacing 



the domains. We believe that the results presented here 
will arise fresh experimental and theoretical interest in 
this problem. 
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